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Abstract 

Cracks are the major vehicle for material failure, and often exhibit rather 
complex dynamics. The laws that govern their motion have remained an 
object of constant study for nearly a century. The simplest kind of dynamic 
crack is a single crack that moves along a straight line. We first briefly review 
current understanding of this "simple" object. We then critically examine the 
assumptions of the classic, scale-free, theory of dynamic fracture, and note 
when it works and how it may fail if certain of these assumptions are relaxed. 
A number of examples is provided, where the introduction of physical scales 
into this scale-free theory profoundly affects both a crack's structure and the 
resulting dynamics. 



1 Introduction 

1.1 Physics of Cracks 

The ability of solids to withstand mechanical forces is one of their fundamental 
properties. Solids divide roughly into two classes, brittle and ductile. The brittle 
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solids break easily, catastrophically. Ductile solids flow or bend before they break 
and are more resilient. There is a simple thermodynamic definition of an ideally 
brittle solid. 

Consider any solid sample of height b whose top and bottom boundaries are 
displaced by distances ±5/2, as shown in Fig. Q] If the total displacement 8 is 
sufficiently small, then the restoring force F on the boundaries will be F = CA8/b, 
where C is an elastic constant, and the change in energy of the solid must be 

AE = ^CA8 2 /b. (1) 

For all solids one can also define the surface energy 7, where 2y is the minimum 
energy per area needed to separate bonds and cut the solid into two pieces. For a 
given solid 7 may depend upon the angle of the cut, but we will assume for simplic- 
ity that the solid is oriented so that 7 has a minimum along the plane perpendicular 
to F. Then there is a critical displacement 8c for which the solid will have lower 
energy by separating into two pieces than by continuing to withstand tension. This 
critical displacement is given by 

2 7 A = \cA8ljb (2) 



^8 G = ^Ayb/C ^F/A = ^AyC/b. 

and is called the Griffith point. In the limit of large systems, where b — > °°, the 
force per area goes to zero. This means that the approximation of Eq. Q] is as 
accurate as desired. Thus, there exists a universal relation between energy, force, 
and boundary displacement shown in Fig. [H and given by 



\CA8 1 lb for 8 < 8 G _ (CA8 /b for 8 < 8 G 
\CA8ljb for5><5 G . ~ [O for 8 > 8 G 



For large enough samples this relation is exact. 

A solid that actually obeys Eq. [3] is ideally brittle. Ceramics such as silicon, 
and some glasses come close. However, for most materials the whole picture is 
misleading or simply wrong. It is not enough for thermodynamics to say that a 
solid can lower its energy by dividing into two pieces. There must actually exist 
some physically feasible way for the separation to take place. Accomplishing this 
separation is the role of cracks. Cracks are nonequilibrium propagating dissipative 
structures. Initiating at a weak spot, they travel across solids under mechanical 
stress and separate them into multiple pieces. The questions of which solids are 
brittle, which are ductile, and how much force or energy truly is needed in order 
to cause a solid to fail come down to the questions of when cracks propagate and 
how. 
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1.2 Scope of this article 

Two of us wrote a review 10 years ago that laid out the basic mechanical theory of 
cracks and then described conditions under which the tips became unstable [ 1]. The 
opportunity to write a second review raised a question of how to cover additional 
topics in an article that would still be fairly self-contained. We decided to focus on 
subjects that on the one hand represent recent developments, but on the other hand 
are in some ways more elementary than those we discussed before. We examine 
both experimental and theoretical issues that arise when one checks carefully if the 
theory of fracture is complete. 

The fundamental theory of fracture, called Linear Elastic Fracture Mechanics 
(LEFM), is now a well established body of knowledge. The mathematical foun- 
dations were carefully laid by investigators such as Irwin [2], Rice [3] and Willis 

and are documented in textbooks such as those by Broberg Q, Freund O, 
and Slepyan Q. By itself, this theory of fracture does not have a length scale 
within it, and it must be supplemented with some physical information about what 
makes cracks start to propagate that can either be obtained from experiments, or 
deduced from additional calculations. Here we will discuss a number of different 
cases where scale-free fractures interact with phenomena involving specific physi- 
cal scales. 

In contrast with physical theories for solid properties such as electrical con- 
duction, physical theories of fracture are still in a primitive state. The ability to 
predict material properties based upon knowledge of atomic constituents or other 
underlying properties, is limited. We have assembled recent advances in describing 
detailed processes near crack tips in the hopes that, eventually, the understanding 
of crack tip mechanics will reach the point where crack direction and speed come 
under reliable experimental control. It may then be possible to design materials so 
that cracks begin to move at designated thresholds. At the present time these goals 
are not in reach, as we describe in our concluding section on conclusions and open 
questions. 

2 Fracture Mechanics 
Linear Elastic Fracture Mechanics 

When a crack exists in a material, any externally applied stress undergoes a very 
large amplification at its tip. This was first noted by Inglis [8], who demonstrated 
amplification of the stresses at the tip of an elliptical hole in an otherwise uniformly 
stressed sheet. Irwin and Orowan JUEl later showed that if a crack exists in a sheet 
assumed to obey linear elasticity, the stress field Oy actually becomes singular at 
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Figure 1 : (A) Universal relation between energy, force and displacement for ide- 
ally brittle solids in tension. The energy reaches a plateau and force vanishes at 
the Griffith point do- (B) Geometry of sample for which these universal relations 
between energy, force, and displacement can be observed. The top and bottom of 
the sample are gripped rigidly and moved apart by distance 8/2. A crack through 
the middle of the sample provides a way for it to break in half. (C) Closeup view 
of universal features of displacements and stresses near the tip of the crack in (B). 
Very close to the crack tip linear elasticity breaks down; sometimes it can be re- 
placed by other theories described later in this paper. 

its (infinitely sharp) tip. The stress (see Fig. \VP) at a distance r from the tip takes 
the universal singular form Of/ oe K/y/r where the coefficient K of the singularity is 
called the stress intensity factor. Depending on the symmetry of loading, there can 
be three independent fracture modes 0121 El, Mode I - fracture in pure tension, 
Mode II - fracture in pure in-plane shear, and Mode III - fracture by out-of -plane or 
"tearing". The stress field due to each of these modes has the same universal form, 
and there may be, in general, three independent stress intensity factors. Here, we 
will only consider fracture under pure tension (Mode I). 

The existence of a singular stress is, itself, insufficient to make a crack prop- 
agate, since the singularity is always cut off at least at the atomic scale, where 
it may or may not be large enough to snap bonds. A better way to determine if 
cracks can move is to note that new surface must be created and this requires en- 
ergy. Defining the energy required to create a unit area as the "Fracture Energy", T, 
a simple generalization of Griffith's idea [11] suggests that cracks begin to move 
when the potential elastic energy per unit area released by a unit extension of a 
crack becomes equal to Y (cf. Eq. O. Even if this condition of energy balance 
predicts the onset of motion, once a crack starts to move the kinetic energy due to 
material motion around the moving crack has to be taken into account. Accounting 
for this energy led to the first equation of motion for a crack, which Mott initially 
obtained by dimensional analysis lPT2l[T3Tl . This equation predicted that in a two- 
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dimensional medium of infinite extent, a crack should continuously accelerate as 
a function of its instantaneous length to a limiting, but finite, asymptotic velocity. 
This speed limit for a crack is due to physical constraints on energy transport. The 
crack tip requires energy to extend. If this energy has to be transported to the tip 
from remote parts of the medium by elastic waves, the speed of a crack must cer- 
tainly be limited by the propagation speed of those waves. This asymptotic limit 
can be reached by either increasing the amount of energy driving the crack or by 
reducing the fracture energy T to zero. Stroh [14] noted that for T — > 0, a crack 
propagating at its asymptotic velocity is equivalent to a disturbance moving along a 
free surface and therefore predicted the crack's limiting velocity to be the Rayleigh 
wave speed, Vr, which is the highest speed that a wave can move along a free 
surface. 

These intuitive ideas were later lfl5l[T6l l6ll4"ll3l shown to be rigorously correct, 
when a quantitative theory of dynamic fracture was developed. This complete 
theory is called Linear Elastic Fracture Mechanics (LEFM). Assuming that the 
medium always obeys linear elasticity, the theory predicts that the general form for 
the stress singularity at the tip of a crack as a function of the distance r and angle 
6 from the tip of a crack moving at velocity v is given by 

where fjj(v,6) is a known universal function and K can be explicitly calculated 
for any externally applied loads. K has the dimensions of stress x length 1 ' 2 , and 
thus must depend on a macroscopic geometric lengthscale. For sufficiently large 
samples, there is always a region around the crack tip where the singular term in 
Eq. @] dominates any other contribution to the stress field. This fact can then be 
used to find the energy release rate, G, defined as the energy per unit extension 
per unit width that is flowing into the tip of a crack. G is related to K by 0: 
G = (1 — v )K A(y)/E, where v and E are, respectively, the Poisson ratio and 
Young's modulus of the medium. A(v) is a universal function of the instantaneous 
crack velocity v. With this expression in hand, the general equation of motion 
for a crack can be obtained by using the energy balance criterion, equating G to 
r. Energy balance, however, needs to be supplemented by an additional condition 
that tells us in what direction a crack will move. In general, no such first-principles 
condition exists. One commonly used assumption, called the "principle of local 
symmetry" ifTTl [T8l . tells us that a crack will locally align itself in a direction 
so as to negate any shear component at its tip. For quasistatic ("slow") fracture 
this idea has been shown to be justified by symmetry considerations (I9l |20| and 
has been used, for example, to quantitatively describe ETl 1221 1251 1241 1251 1261 1271 
l28l I29I experimental measurements of slow oscillating cracks in heated strips QUI 
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[3TTl . Other closely related path selection criteria were proposed and tested, see for 
example ||32l . Although the "principle of local symmetry" is often assumed for 
rapid cracks, there is, however, neither fundamental or experimental justification 
for this (or any other) crack path selection criterion. 

An additional tacit assumption is also used; that all of the dissipation in the 
system occurs within a "small" region surrounding the tip of the crack. This means 
that all of the complex dissipative mechanisms, inherent in the fracture process, 
have to occur in a region for which the surrounding stress field is entirely domi- 
nated by the singular stress field described by Eq. HI This condition, called the 
condition of small-scale yielding, is, in essence, the definition of a brittle mate- 
rial - which extends the concept of an ideal brittle material (cf. Fig. [IK). The 
small zone surrounding the tip in which all nonlinear and dissipative processes oc- 
cur is called the process zone. Thus, T need not necessarily be the energy cost to 
break a single plane of molecular bonds - it embodies all of the complex nonlinear 
processes that are driven by the putatively singular stress field near the crack tip. 
Small-scale yielding is one of the foundations of fracture mechanics, as it ensures 
the equivalence of the global energy balance criterion used by Mott with the locally 
formulated equation of motion: 

G = r=(l-v 2 )K 2 (l,v)A{v)/E. (5) 

Freund [ 6 ] demonstrated that for straight semi-infinite cracks in infinite plates, un- 
der a variety of loading conditions 0, the dynamic stress intensity factor K(l,v) 
has a separable form K(l,v) = K(l)k(v), where K(l), which can be explicitly cal- 
culated, is solely dependent on the external loading configuration and the instan- 
taneous value of the crack length, /. Both k(y) and A(v) are universal functions 
of only v. Freund |6] showed that, to a good approximation k(v)A(v) wl - v/Vr 
yielding 

r_G«(l^Q.(l-v/ft) (6) 

t, 

Inverting Eq. © yields the following prediction for an equation of motion: 

YE 



Vr 



1 



(7) 



(1 

which is identical to the equation obtained earlier by dimensional analysis for a 
crack in an infinite plate. 

Does this equation of motion work? Sharon et al. |[33Tl conducted a quantita- 
tive test of Eq. [7J where the detailed dynamics of rapidly propagating cracks were 
measured in both PMMA (Plexiglas) and soda-lime glass, two brittle amorphous 
materials. The experiments measured the instantaneous values of / and v under 
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loading conditions whose geometry enabled a precise calculation of K(l). The val- 
ues of r were derived using Eq. |6]for widely different experimental conditions. T 
is a material property, but need not be velocity independent and, indeed, is gener- 
ally a function of v. If Eq. |7]is valid, the values of T(v) derived from Eq. |6]using 
measurements of / and v, should all collapse onto a single curve for each material. 
The results of these experiments were surprising. As we see in Fig. |2j the ex- 
periments for both materials indeed collapse onto single curves for both glass and 
PMMA at crack velocities that are less than about OAVr in each material. Thus, 
Eq. |7] is indeed in excellent quantitative agreement for v < OAVr for both glass 
and PMMA, two materials whose micro-structure is entirely different. These re- 
sults agree with earlier results of Bergqvist 041 and Ravi-Chandar et al (35 ] that 
showed that cracks appear to behave in a way that is consistent with LEFM, for 
relatively low crack velocities. 

What happens when v > 0.4Vr? At a critical velocity, v c ~ OAVr a single 
crack state become unstable 1381 [Q [39l to one in which the main crack is continu- 
ously sprouting frustrated microscopic branched cracks, as shown in Fig. [3t. This 
"micro-branching" instability causes rapid oscillations in the instantaneous crack 
velocity and the formation of non-trivial structure on the fracture surface. There 
is currently no first-principles theory for the origin of this and other instabilities, 
although many serious attempts have been made rt40ll4Tl l43l l43l l48l 
@9l|50l[5ll|52l|53l[5ll|55l[56l[53 It seems clear, 

however, that the key does not lie within the framework of LEFM. In fact, once a 
simple single-crack state is lost, we have little fundamental understanding of frac- 
ture dynamics. Some examples of the complexities that such non-trivial states can 
produce are shown in Fig. [3] 

Does LEFM break down for v > OAVr if instabilities are suppressed? Recent 
experiments Ir36ll67ll37i , 68 ] show us that, as long as its underlying assumptions are 
met, LEFM is entirely accurate. For example, the displacement fields surrounding 
a crack's tip can be obtained by integrating Eq. H] In particular, displacements 
normal to the crack's propagation direction (6 = 0), predicts a parabolic crack tip 
opening u y (r < 0, 6 = 7t) °< K(l,v) • yfr behind the crack (Fig. [QC-top), whose 
curvature is determined by K(l,v). Crack instabilities have been suppressed to up 
to v ~ O.Wr in soft brittle gels (671 . In these materials, Vr is 1000 times lower than 
in standard materials. This enabled direct measurements of the crack tip profile 
of rapidly propagating cracks IT371 . These experiments show that, as long as we 
are not too close to the crack tip, LEFM provides excellent quantitative agreement 
with measurements. Values of T(v) collapse onto a well-defined function over a 
very large range of velocities, as Fig. [2£ shows, when calculated using values of 
K(l,v) obtained by measuring the crack tip curvature. 
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Figure 2: Fracture energy, T, derived using Eq. [6j as a function of the mean ve- 
locity for soda-lime glass (A) and PMMA (B). Top: Velocities of cracks driven 
with different initial stresses and initial crack lengths as a function of their instan- 
taneous lengths. The data were smoothed over a 1mm length, to filter out velocity 
fluctuations. The dashed lines indicate the highest values of v used to derive the 
corresponding T(v) (lower plots). These are 350 m/s (0.38 Vr) in PMMA and 1325 
m/s (0.40 Vr) in glass. The observed data collapse to a single (material-dependent) 
function, T(v) is a quantitative validation of the equation of motion derived for 
a single crack. Data were taken from 11331 . Values of Vr for PMMA and soda- 
lime glass are, respectively, 930m/s and 3,300m/s. (C) T(v), derived by using the 
parabolic profile of the crack tip predicted by LEFM. These measurements were 
made possible by lt3"6ll37l using soft brittle gels, where Vr = 5.3m/s. 
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2.1 LEFM - when does it work? 



Sometimes LEFM provides an excellent quantitative description of both a crack's 
motion and of the elastic fields surrounding its tip. This happens whenever the, 
at times subtle, fundamental assumptions underlying LEFM are valid. Let us list 
these assumptions: 

• Only the dynamic behavior of straight single-crack states are accounted for. 
Once more complex dynamic states (e.g. micro-branching) occur, LEFM 
fails to accurately predict their behavior. In fact, there is currently no first- 
principles criterion for what path a crack must take. 

• LEFM assumes that linear elasticity is valid away from the crack tip. In 
principle, for large enough samples, this assumption can always be made 
as accurate as desired, but in practice for samples of fixed size it must be 
checked. 

• Small-scale yielding is observed 44> all nonlinear and/or dissipative processes 
are assumed to occur in a region of negligible size at a crack's tip. As LEFM 
is a scale-free theory, this tacitly means that the scales at which these com- 
plex processes occur are negligible. We will see that this tacit assumption is 
not necessarily justified; the existence of additional scales within the crack 
tip region can and does have a large effect on crack dynamics. 

• A tacit assumption concerning "energy balance" is that energy has to flow 
into a crack's tip, otherwise a crack cannot propagate. 

In the current "state of the art" in fracture we have little fundamental understanding 
of the complexities of a crack's motion once a crack decides to stray from a straight 
path or become unstable. These questions are hard ones, and although there is a 
number of possible directions to take, we have no clear answers. In the following 
sections we have, instead, chosen to describe what we currently know about "sim- 
ple" straight cracks. We will do this by describing a number of interesting cases 
where theory and experiment move beyond limits constraining LEFM. By doing 
this we hope to provide both insights into the character of the "simple" crack state, 
and to, perhaps, provide a springboard for understanding more complex fracture 
dynamics. 
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Figure 3: (A) A cross-hatching instability in slow hydrogel Mode I fracture due to 
straight macrosteps on the crack surfaces l69l . (B) A high-speed crack tip oscil- 
latory instability in Mode I fracture of an elastomer gel l67l . (C) A side view of 
the side-branching instability in rapid Mode I fracture in Plexiglas (left column) 
and the accompanying crack surfaces (right column) iTTOl [Tl. (D) Complex fea- 
tures on the fracture surfaces of a continuously twisting mixed Mode I+III crack in 
Plexiglas EQ. 

3 Tales of Scales 

3.1 The Importance of Small Scales near a Crack's Tip 

LEFM is a scale-free theory. Scales, however, must be introduced (either tacitly 
or explicitly) to obtain a closed theoretical description of a crack's dynamics. A 
necessary scale, for example, is inherent in the idea of the process zone, which 
both cuts off the singular fields given by Eq. [4] and provides a scale below which 
dissipation and material nonlinearity take place. A central question in fracture is 
whether the details of these small scales are important. Is a proper description of 
these scales needed to describe how a crack behaves, or do the universal properties 
of fracture indicate that the details of the fracture processes that take place at these 
scales are irrelevant and simply serve to determine a material's fracture energy. 
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There is a number of hints that the character of these small scales plays a major 
role in fracture dynamics. A classic way to regularize the continuum singularity 
at the tips of cracks is to introduce cohesive zones, fictitious forces pulling the 
crack's faces together behind the tip in just the right way so that the singularity 
of the stress field disappears [72] [73]]. In cohesive zone models, these fictitious 
forces acquire equations of motion so that they can travel along with the crack, 
cancelling out the singularity as they go. Results of working out such models are 
disquieting. It is difficult to obtain stable straight-ahead propagation of a single 
crack at any velocity [4911741 . As noted by Lobkovsky and Langer [74], this entire 
class of models "seem[s] to be highly sensitive to details that ought to be physically 
unimportant". 

This point of view is reinforced by numerical simulations. For example, fi- 
nite element calculations of dynamic cracks 11731 with cohesive zones can provide 
quantitative agreement with experimental descriptions of both stable cracks and 
cracks undergoing the micro-branching instability described in Fig. However, 
it turned out that even the existence of micro-branching in numerics depends on 
the mesh size used in the simulations, and the experimental phenomena disappear 
in the limit where the mesh size goes to zero IT761 . 

In this section, we will first describe three approaches in which a well-defined 
length scale at the crack's tip occurs naturally. Two of these are continuum theories 
in which scales are formed either by (1) relaxing the assumption that materials 
are elastic up until the point where dissipation occurs or (2) the introduction of 
a phenomenological scalar field that represents the material damage that occurs 
around the tip. Both of these approaches define a new dynamic scale at the tip. 
A third class of models will be described in which the medium is an ideal crystal, 
which therefore possesses an intrinsic length-scale, the lattice size. We will see 
that in all of these approaches the existence of a length scale at the crack tip has 
important consequences on both the dynamics and form of propagating cracks. 

We then show that an entire new class of single-crack solutions can arise when 
the Griffith condition is supplanted by the introduction of new energy scale. These 
solutions, which are wholly physical, surpass the speed limit set by LEFM by vio- 
lating its assumption about how energy is conserved near a crack tip. 

3.2 Small Scales Dominated by Nonlinear Elasticity 

An essential feature of crack propagation is the high concentration of deformation 
in the immediate vicinity of a crack's tip, as given by LEFM's y^r-singular fields. 
LEFM, however, is inherently confined to small material deformations. Therefore, 
any attempt to understand the physics of the process zone involves first the question 
of how and where LEFM breaks down near the tip of a crack. To address this 



11 



question, consider the interatomic potential by which particles interact. Adopt now 
the simplified assumption that the material's response until the point of failure is 
completely contained within this interaction potential (i.e. neglect features such as 
dislocation motion, or changes in the underlying geometry of amorphous solids). 
Linear elasticity is described by the harmonic approximation about the equilibrium 
state at the bottom of the potential well. Plastic deformation, damage evolution 
and finally fracture correspond to interatomic separations far from the minimum 
of the potential, where the interaction energy decreases significantly. It is both 
natural and physically sound to assume that in numerous materials before one of 
these irreversible processes occurs, the harmonic approximation will first break 
down via reversible nonlinear deformation. Therefore, the first physical process 
that intervenes when LEFM breaks down should result from corrections that arise 
from nonlinear elasticity. The first generic nonlinear elastic contribution would be 
expected to come from quadratic corrections to linear elasticity. 

To model this behavior at a macroscopic level, we expand a general stress vs. 
displacement-gradient relation up to second order 

O~lldu-fi(du) 2 + 0l(du) 3 }. (8) 

For simplicity, we suppress the tensorial nature of the quantities involved and 
let a and du schematically denote stress and displacement-gradient (strain) mea- 
sures, respectively (u denotes a displacement and d schematically denotes a spatial 
derivative); /I and jl denote the first and second order elastic moduli, respectively. 
We term the theory that is based on this expansion "weakly nonlinear fracture me- 
chanics" ll36l 1771 l68l . as higher order nonlinearities are neglected in Eq. [8] To 
explore the implications of Eq. [8] on fracture dynamics, we schematically write the 
momentum balance equation as 

do = , (9) 

where, for simplicity, we omit the inertial term on the right hand side. We then 
expand the displacement u in powers of the magnitude of the (small) displacement- 
gradient £ ~ \du\ and to second order obtain 

M ~ £S (i) + £ 2~(2) _ M (i) +M (2) _ (10) 

Substituting Eq. [T0]in Eq. [9] and expanding to first order, we obtain 

Hd 2 u {l) =0, (11) 

which is our schematic analog of the Lame equation f6 j. LEFM tells us that this 
equation, together with the traction-free boundary conditions on the crack faces, 
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leads to the following asymptotic behavior near the tip of a crack 1161 



d M (l)^JL (12) 

where K is the stress intensity factor and r is measured from the crack tip, cf. Eq. 
IU In terms of the displacement-gradients Eq. |9j to second order, becomes 

lid 2 u ( V-fid(du^) 2 = 0. (13) 

This is analogous to a Lame equation for i/ 2 ' with an effective body-force which 
scales as 

jG a(a M (i))2^^!. ( i4) 

/I T 

By inspection, Eq. (fT3l admits a solution of the form 



du^ ~ q_ . ( i5) 



We therefore see that simply accounting for quadratic contributions to material 
stress-strain relations already leads to the interesting conclusion that the first non- 
linear correction to the asymptotic LEFM fields is characterized by a 1 jr displacement- 
gradient singularity. The simple scaling-like considerations described above cap- 
ture the essence of the solution, although the complete dynamic (i.e. with inertial 
effects) weakly nonlinear solution is, of course, much more mathematically in- 
volved 061 1681 . For example, the complete solution, which depends on the tip 
polar coordinates r and 6, shows that 1 Jr displacement-gradients terms arise from 
both r-independent and log (r) displacement contributions; the latter has an im- 
portant implication as it yields significant distortion of the parabolic crack profile 
when the tip is approached. 

Having addressed the question of how LEFM breaks down under the stated 
assumptions, we proceed now to consider the question of where this happens in 
space. According to the preceding discussion, LEFM breaks down when the sec- 
ond order correction du^ becomes non-negligible with respect to du^\ A plausi- 
ble estimate for this condition is obtained from the condition that the second order 
nonlinear contribution to the stress-strain relation becomes larger than 10% of the 
leading, linear, term. Thus, 

M»^ 0A - (16) 

This immediately implies that higher order corrections like 1 / r 3 / 2 are negligibly 
small in this region. As the strains surrounding any crack tip are approaching a 
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Figure 4: (A) Left: An experimental picture of a crack, with a parabolic LEFM fit, 
cf. Fig. [2] Middle: A zoom in on the crack tip. Right: A comparison of the weakly 
nonlinear theory (solid line) with the measured crack tip opening profile (circles). 
The LEFM parabolic fit is added (dashed line), (x' ,y') are the coordinates in the 
deformed (laboratory) frame. (A) A comparison of the weakly nonlinear theory 
(solid line) with the measured strain £ yy along the symmetry line ahead of the crack 
tip, for three crack velocities. The LEFM predictions are added (dashed lines). See 
text for more details. 



mathematical singularity, this condition must always be satisfied at some scale. 
Denoting the scale at which Eq. [16] is first satisfied by i„i and using Eqs. [12] and 
[T5l we obtain 

This simple analysis reveals the emergence of a new length-scale associated with 
weak elastic nonlinearities. It is important to note that l n i is dynamic as its value 
increases with the energy driving the crack ~ K 2 / /I or, equivalently, with a crack's 
velocity. 

Is the scale i n i in Eq. [I6]physically relevant? Nonlinear elasticity in large-scale 
measurements is often not readily apparent. Brittle materials nearly always appear 
linear elastic until the onset of irreversible deformation. These measurements can 
be, however, quite deceiving. Failure in these tests is not due to intrinsic material 
behavior, but generally takes place via the activation of cracks or defects that are 
already in place within the material. On the other hand, cracks generally propagate 
into locally virgin material. The deformations near the tip of a crack are large 
enough to be well within the nonlinear elastic regime. As defects are ubiquitous 
in most macroscopic materials, these effects may be nearly undetectable in large- 
scale measurements. This does not mean, however, that they are either unimportant 
or do not commonly occur. 



14 



With this in mind, irreversible (plastic) deformation may commonly be pre- 
ceded by elastic nonlinearities as the displacement-gradients increase. Thus, the 
plastic deformation scale £ p i is expected to satisfy 

tnl>tpl. (18) 

For strongly nonlinear elastic materials (gels, rubber etc.) we expect a scale sepa- 
ration of the form l n i^>t p \, while for other materials we cannot rule out the pos- 
sibility that l n \ is not very different from £ p i. In the former case, the small scales 
near a crack's tip are dominated by nonlinear elasticity, while in the latter case, 
dissipation may play a more decisive role. 

The complete solution of the weakly nonlinear equations, which generalizes the 
simple considerations above, was obtained in Refs. Il36ll68l . Can these solutions 
be observed in experiments? Generally, this is a near-impossible task since cracks 
propagate at near sound-speed velocities with (at best) /xm-scale crack tip open- 
ings. In addition, the second order elastic constants are generally not well known, 
since in macroscopic stress-strain tests of brittle materials failure, due to the propa- 
gation of intrinsic cracks, occurs well before sufficient strains are reached. Recent 
experiments, however, circumvented these problems by utilizing brittle elastomers 
(gels) [37]. While the characteristic fracture phenomenology of these soft materi- 
als is identical to that of other brittle amorphous materials [78], their sound speeds 
are 1000 times lower. Furthermore, as these materials are very compliant, but 
tough, they undergo large strains prior to fracture. This allows precise measure- 
ments of their second-order elastic constants. These properties enabled the first 
direct experimental measurements of the deformation near the tip of rapid cracks. 
A summary of the comparison of these measurements with the explicit predictions 
of the weakly nonlinear theory is presented in Fig. [4] On the left hand side of the 
upper panel the LEFM parabolic crack tip profile is shown, cf. the discussion of 
this point in Sect. [2] and Fig. [2] In the middle of the upper panel a zoom in on 
the crack tip region is shown, where marked deviations from the LEFM parabolic 
profile are observed. In 11371 it was shown explicitly that this deviation is due to 
elastic deformation and not due to an irreversible one. 

On the right hand side of the upper panel we show the predictions of the com- 
plete weakly nonlinear solution for the tip profile. We observe that it agrees well 
with the LEFM parabolic form at large scales (i.e. small deformation), but both 
qualitatively and quantitatively (with no adjustable parameters, since all elastic 
properties can be measured in this material) captures the deviation from it on 
smaller scales (below ~ 200/xm in this example). The deviation is a direct con- 
sequence of the log (r) contribution to u x and directly demonstrates the existence 
of the new terms predicted by the weakly nonlinear theory. 
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In the lower panel of Fig. @]we present a comparison between the weakly non- 
linear theory and the measured deformation ahead of the crack tip. We focus on 
E yy = dyUy along the symmetry line (6=0) ahead of the crack for three propaga- 
tion velocities, including one that approaches the limiting speed, Vr. For the two 
lowest velocities we observe marked deviations from the 1 / y/r fields of LEFM, 
deviations that are accurately captured by the weakly nonlinear theory (cf. 061 for 
more details). These results again explicitly demonstrate the validity of the weakly 
nonlinear theory with its l/r singularity. Moreover, the scale l n i in which LEFM 
breaks down via elastic nonlinearities is captured by the theory. Its value, which is 
in the mm range, can be easily read off from the figure. 

On the right hand side of the left panel we present the same comparison for a 
crack propagating at v=0.78c s . At this very high velocity, higher than second order 
nonlinearities are needed to accurately describe the data. This comparison, how- 
ever, highlights an important point. According to LEFM, a velocity exists (0.73c s 
for an incompressible material) where % changes sign from positive to negative. 
In this range of velocities Eyy predicted by LEFM is small (due to kinematic func- 
tions [6j |37j [36 ] that enter as coefficients in the exact analog of Eq. [121 but do not 
appear in the simple derivation used above). As a result, in this range of velocities 
the second order term in £ yy becomes the dominant contribution. Due to this contri- 
bution, as seen in the figure, no change in sign occurs. Thus, the weakly nonlinear 
theory retains our basic physical intuition about fracture; material points straddling 
the symmetry line must be separated from one another (e yy > 0) to precipitate frac- 
ture. This contrasts with the LEFM prediction in this range of velocities where 
e yy <0 (cf. dashed line in the figure). 

In summary, the weakly nonlinear theory of dynamic fracture describes how 
LEFM breaks down in materials where the scale defined by elastic nonlinearity is 
larger than the dissipative scale. In addition to accurately describing detailed mea- 
surements of the deformation near tips of very rapid cracks 061 1681 . this theory 
may also have further profound implications for understanding crack tip instabili- 
ties H][36l. This theory, by itself, does not account for either near-tip dissipation 
or directly offer crack path predictions. The length-scale £„i for the breakdown of 
LEFM may, however, play a central role in explaining symmetry breaking insta- 
bilities. There are some tantalizing hints that crack tip instabilities may occur at 
approximately the length-scale l n \ that emerges in this theory. For example, l n i at 
high velocities correlates well with the geometry-independent wavelength of crack 
path oscillations discussed in Il67ll66l . This promising line of investigation should 
be further explored. 
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3.3 Phase Field Models 



LEFM accurately predicts the transport of linear elastic energy from the large 
scales, where external forces are applied, to the small scales near the tip of a crack 
(the process zone), where energy is dissipated in fracture. It is assumed that these 
scales are coupled solely through the y/r - singular field. This approach, however, 
is inherently deficient in at least three important aspects: 

• The crack tip velocity v cannot be determined self-consistently within LEFM, 
as energy dissipation within the nonlinear near-tip region is not described in 
this framework. 

• The path selected by a crack's tip remains undetermined in the framework of 
LEFM unless supplemented by a path selection criterion (e.g. the principle 
of local symmetry). Hence crack tip instabilities cannot be explained by 
solely LEFM. 

• LEFM is a scale-free theory and hence any phenomenon that involves a non- 
geometric length-scale is beyond its scope (e.g. some crack tip instabilities). 

All of these difficulties are directly related to the need to account for the physics 
at the small scales near the tip of a crack, where LEFM breaks down. A compre- 
hensive theoretical framework should explain how LEFM breaks down near the tip 
of a crack, predict crack tip shape and velocity selection, quantify energy dissipa- 
tion near the tip of a crack and account for the path selected by cracks. Crack tip 
instabilities are expected to naturally emerge from such a theory. 

A class of phenomenological phase field models has recently been developed to 
approach this problem from a continuum perspective f79l 180118 Ill82ll83ll84ll§5ll86l 
|29l[87l. In this approach, the fracture of materials is described as the gradual accu- 
mulation of damage near the tip of a crack. Damage is quantified mathematically 
by a scalar "order parameter" or "phase field" (j>(x,t), where = 1 corresponds to 
an intact material and = to a fully broken one. The physical interpretation of 
the phase field as quantifying local damage is, in fact, not necessary and most 
probably will not be advocated by all practitioners of this approach. The phase 
field can be viewed as a mathematical interpolation between LEFM and material 
failure, where = 1 is favored when the linear elastic strain energy density S e i is 
smaller than some threshold S c , <§ e \ < S c , while = is favored when S>S C . A 
similar approach proved to be very fruitful in the context of non-equilibrium crystal 
growth, when used to predict solidification micro-structures. 

The mathematical formulation of these ideas starts with the following expres- 
sion for the potential energy density as a function of and the displacement field 
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u 

^ = f |V0| 2 +g(0)(4/-^ c ) • (19) 

$ e i is the linear elastic energy density, which is given in terms of linear elastic strain 
tensor s, as 

<£z = ^e£ + juej , e i j = ^{dju i + d i u j ) . (20) 

A, jU are the Lame constants. Two of the phenomenological quantities introduced 
by phase field models appear here: K, a coefficient with dimensions of energy 
per unit length and g(<j>), a monotonic function that quantifies the degree of strain 
softening due to material degradation. g(<p) satisfies g(0) = 0, g(l) = 1, g'(0) = 
g'(l) = 0, and lim0^o^(0) ~ 0", with a > 2. The latter condition was shown to 
ensure the existence of crack-like solutions I8T1 . We note that <p in Eq. [l9]does not 
break the rotational invariance near the crack tip, unlike lattice or cohesive zone 
approaches that inevitably introduce some near-tip anisotropy. This anisotropy 
is questionable when amorphous materials, which account for a large part of the 
available experimental data on rapid fracture 12, are considered. Note also that 
Eq. [19] is symmetric with respect to tension and compression, though it is obvious 
that fracture occurs locally under tension. This issue was addressed in ll83l by 
introducing an ad-hoc asymmetry between tension and compression. 

The equations of motion for a phase field model are variationally derived from 
the total potential energy E = j SdV and read ll82ll83Tl 

pd ttUi = d ] ^--— i Hij^djUt, (21) 

1 M = dj^--^, Zj^djt, (22) 

where p is the mass density. Translational invariance implies that $ is independent 
of u, i.e. that |^ = . Noting that j^- = a ( y, Eq. [21] is, in fact, the usual linear 
momentum balance equation, pd tt iii = djOij. 

The dynamics of <p in Eq. [22]are assumed to be dissipative and hence first order 
in time. The last phenomenological parameter %, needed to close the theory, has 
the dimensions of volume per unit energy per unit time. Equations (l2TI) - (|22l are 
four coupled nonlinear equations that can be solved numerically when the external 
boundary conditions are specified. The simple choice ofg(0)=40 3 — 30 4 , which 
satisfies all the required conditions mentioned above, is commonly used in the lit- 
erature ||8T1 |82l [83l . It is important to note that the usual traction-free boundary 
conditions on the crack faces @ are not needed here since these are automatically 
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Figure 5: (A) Branching patterns in a two-dimensional out-of -plane (Mode III) 
phase field model simulation Il82l . The external loading level is increased from top 
to bottom. (B) Fast sinusoidal crack tip oscillations in a two-dimensional in-plane 
(Mode I) phase field model simulation under bi-axial loading ll83l . (C) Quasistatic 
nonlinear crack tip oscillations in a two-dimensional in-plane (Mode I) phase field 
model simulation under thermal loading [29]. (D) Complex crack front morphol- 
ogy where crack breaks into segments in a three-dimensional phase field model 
simulation under mixed mode loading (in and out-of-plane, Modes I+III), courtesy 
of A. Karma and A. J. Pons. 



satisfied due to strain softening when 0^0. Moreover, fracture in the phase field 
models is described by a diffuse interface, which is a significant technical advan- 
tage, since one does not have to track a moving boundary, which is computationally 
expensive. The crack faces in the phase field models are rather arbitrarily defined 
as a contour of fixed (j), usually = 1/2 lf8Tl 1821 l83t 

A simple dimensional analysis, based on the phenomenologically introduced 
parameters K, S c and %, shows that the size of the region where LEFM breaks 
down, the process zone, scales as 




(23) 



the fracture energy T scales as 



r-V^, (24) 

and the timescale for energy dissipation near the tip of the crack scales as 

T ~ — . (25) 
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Therefore, in light of the discussion at the beginning of this Section, the phase 
field models include a regularizing cut-off length-scale £ and near-tip dissipation, 
resulting in a well-defined and self-consistent mathematical description of a frac- 
ture problem. Thus, crack tip shape, path and velocity selection are expected to 
emerge naturally from the numerical solution of Eqs. [2111221 in this class of mod- 
els. We note that the variational structure of the phase field models allows a natural 
generalization of common energy-momentum tensors and configurational forces in 
fracture mechanics ll87l . 

Several groups have explored numerical solutions of the phase field equations 
in various configurations. The results exhibit a rather rich phenomenology that is 
at least in qualitative agreement with experimentally observed phenomena; these 
include two-dimensional steady-state crack motion in a strip geometry above the 
Griffith threshold R811I85I . branching instabilities (80l|82l[86l, rapid [83l and quasi- 
static [29 ] oscillations and complex three-dimensional crack front dynamics. In Fig 
. [5]we summarize some of these results. 

The phase field approach seems a very good mathematical method for solv- 
ing quasi-static fracture problems that are completely controlled by LEFM and by 
geometric length-scales, cf. QUI . In these cases it provides a proper crack tip 
regularization and dissipation, and offers a computationally efficient way to track 
complex crack configurations. For example, the morphology of a quasi-static crack 
propagating under thermal stresses was satisfactorily described in the framework 
of a phase field model [29]. In this example, all the pertinent length-scales were 
geometric in nature and the crack path was accurately described by the principle 
of local symmetry. On the other hand, in situations where the crack tip physics 
has explicit macroscopic manifestations lf88l l67l l66ll . a number of important open 
questions remains. The phase field methodology does not offer a physically real- 
istic description of the deformation and dissipative processes near crack tips. It is 
therefore difficult to directly relate the phase field <j>(x,t) itself or the phenomeno- 
logical parameters K and % to experimental data. Therefore, this approach cannot, 
at present, offer testable predictions in these situations. By coupling directly to 
linear elasticity, it also does not incorporate elastic nonlinearities (hence does not 
conform with Eq. [T8T >. Future work is needed to further elucidate the relation 
between the phase field models and the physics of fracture. For example, incor- 
porating nonlinear elastic effects as discussed in Sect. 13.21 may be an interesting 
future line of investigation. 

3.4 Atomic Resolution of Small Scales 

In some materials, details of material structure down to the atomic scale play a 
decisive role in determining how details of the fracture process work out. This 
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statement is most evidently true in crystals, which have natural cleavage planes. 
No matter how big a sample of mica may be, it cracks easily along one easy plane, 
and resists fracture along other directions [89]. One way to characterize the effects 
of the atomic scale is to say that T depends upon orientation relative to crystal 
planes, with cusp-like minima in special directions. The surface energy of crystals 
depends upon orientation in precisely this way [90], and fracture energy must have 
the same general behavior. The fracture of crystals can be strongly anisotropic, 
even in a material that at large scales has completely isotropic elastic behavior 
El. 

While fracture energy has the same qualitative behavior as surface energy in a 
crystal, the two are not equal. Fracture energy cannot be less than surface energy 
because by definition surface energy describes the minimum free energy needed to 
induce a solid to separate along some plane. However it can be greater and in gen- 
eral it is. When cracks travel through a crystal, the motion through the periodic unit 
cells excites phonons. In general, any object moving at velocity v through a crys- 
tal with phonon dispersion relation co(k) excites phonons that obey the Cherenkov 
condition [92] co{k) = vk, and cracks are no exception. While this fact has been 
confirmed repeatedly in theoretical and numerical investigations, we are unaware 
of any direct experimental evidence. 

The obvious way to evade exciting phonons so that fracture and surface energy 
coincide is for the crack to move very slowly. Considering this possibility raises 
a number of interesting problems. Thomson, Hsieh, and Rana [93] showed that if 
one takes a crystal at zero temperature and very slowly pulls its edges apart, when 
the stored energy reaches the Griffith point, the crystal does not break. The strained 
crystal is meta-stable. To see why, imagine that a crack moves along a plane in a 
strip loaded as in Fig. [TJ3. Let give the minimum-energy atomic configuration 
subject to the constraint that the crack tip be at position /. The crack location is 
easy to define as a continuously varying function of atomic positions. Finding 
the function <f (/) in practice is demanding l94l 1931 . but in principle one can see 
immediately what shape it has to have, as shown in Fig. [6) First £ (/) must vary 
linearly with /, with positive slope when 8 is below the Griffith point and negative 
slope when it is stretched above the Griffith point. In addition, $ (/) must have a 
periodic component corresponding to different locations in the unit cell. Right at 
the Griffith point the crack tip sees a corrugated potential, is necessarily trapped 
in a local minimum, and does not propagate. The crystal must be stretched by an 
extra amount so that the slope of <§(V) is negative everywhere. Thomson, Hsieh, 
and Rana [93] called this phenomenon lattice trapping. 

This discussion of linear stability does not immediately rule out the possibility 
that a crack could move very slowly just above the Griffith threshold if it could only 
get started. Exact solutions for cracks running in lattices found by Slepyan ll96l 13 
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Figure 6: (A) Schematic view of how the energy of a crystal varies as a function 
of the location / of a crack tip within a fracture plane. The potential energy surface 
is corrugated, and therefore the crack cannot travel spontaneously at the Griffith 
point when the elastic energy recovered per length equals the surface energy cost. 
(B) Exact solution for crack velocity measured relative to wave speed Vr versus 
loading 8 over Griffith loading 8q in a two-dimensional crystal at zero temperature, 
showing lattice trapping, velocity gap, and the point where steady states become 
dynamically unstable. Past the point of instability steady states are technically 
impossible, but the deviations are on an atomic scale that would be hard to detect 
experimentally. 



and Marder [44] show that this cannot happen either. Every time a crack moves to 
a new unit cell it slides down the corrugated potential and causes the atoms there 
to vibrate at frequencies on the order of c/a where c is a sound speed and a is a 
lattice spacing. If the atom vibrates many times it radiates energy which is lost to 
crack motion. Therefore the most energy-efficient crack motion requires the crack 
to get to a new cell during a time on the order of a vibrational period. This means 
that its velocity must be a fraction of the sound speed, and if the velocity drops 
below a threshold it cannot move at all, as shown in Fig. [6j3. Thus there is a finite 
velocity gap for cracks moving in crystals at zero temperature. Either cracks move 
at a finite fraction of the sound speed and spend energy on phonons or they do not 
move at all. 

All of these conclusions must be modified at nonzero temperature. If the energy 
of thermal fluctuations is comparable to the size of the corrugations in Fig. [6K, then 
the corrugations should be washed out and the crack should behave as one would 
expect if the material is continuous down to the finest scales. 

Attempts to follow through on all these ideas and compare with experiment 
has proceeded the furthest in single-crystal silicon, where the situation is not com- 
pletely settled. The minimum energy plane in silicon is (111), and theoretical 
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Figure 7: Comparison of experiment and theory for the fracture of silicon along the 
(111) plane, showing the energy T(v) needed for propagation at various speeds. 
The experiments were carried out at room temperature in silicon single crystals 
11971 . Deegan established an energy gradient in a strip so that the curve could be 
mapped out in a single run. The three upper theoretical curves are molecular dy- 
namics simulations using the Modified Embedded Atom Method (MEAM) at tem- 
peratures of and 77K, and the Inadvertently Modified Stillinger Weber (IMSW) 
potential [98]. Lattice trapping is so large for these potentials that they disagree 
substantially with experiment. However a tight-binding computation of Bernstein 
and Hess is in agreement |[99l . 
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estimates for the surface energy 2/ range from 1 to 2.5 J/m 2 II100I1 . Experimen- 
tal measurements of the surface energy are in the range of 2-5 J/m 2 11001 , but as 
dynamically propagating cracks have been observed with energies as low as 2.25 
±0.25 J/m 2 |[97l , the lower end of the range is most likely. Both experimental 
and numerical estimates of T(v) are available, as shown in Fig. [7] One might ex- 
pect based upon Fig. [6]\ that the detailed relationship between energy available 
for fracture T and crack speed v might depend rather sensitively upon details of 
the potential energy surface, particularly upon the height of the corrugations, and 
this is the case. The mechanical behavior of silicon at the atomic scale is often 
described by classical empirical potentials 11981 , such as the one due to Stillinger 
and Weber fTOll . or the Modified Embedded Atom Method IT02l . All of these 
empirical potentials appear to have corrugations that are too large, and therefore 
produce substantially more lattice trapping and larger velocity gaps than are seen 
in experiment. The potentials are constrained by considerations of crystalline sym- 
metry, fracture energy, and low-energy properties like sound speed. No one has yet 
managed to guess potentials for silicon that gets all these things right and the frac- 
ture properties as well, although the Modified Embedded Atom Method is so far 
the best. However, an elaborate calculation of Bernstein and Hess |[99l that couples 
tight-binding quantum mechanics near the crack tip to empirical potentials further 
away does seem to capture all of the observed experimental features well. It is a 
bit surprising that such effort must be expended to get the energy at which fracture 
begins right within a factor of 2. 

3.5 Supersonic States: a New Energy Scale 

The preceding sections showed some effects of integrating length scales with LEFM. 
An additional assumption of LEFM (cf. Sect. |2]) was that energy must flow into 
a crack's tip to enable it to propagate. Here we will examine the consequences of 
discarding this assumption. 

The Rayleigh wave speed as the upper limit for a crack was a consequence of 
this assumption. One can ask what limits apply to crack speeds if both the tendency 
to instability is suppressed and the above assumption discarded. Suppressing in- 
stabilities can be arranged by having the fracture move along a weak interface, or 
by finding a material which for some reason does not permit crack branching to 
occur. Cracks travelling along weak interfaces have been of greatest interest in the 
study of earthquakes. Andrews [103] showed numerically that shear cracks on a 
weak interface could travel faster than the shear wave speed. This problem was 
studied further by Burridge, Conn, and Freund II104II . who showed that LEFM per- 
mits shear cracks to travel precisely at a speed given by y/2 times the shear wave 
speed, but that the range of possible speeds is broadened by cohesive forces near 
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the crack tip. 

Since we have maintained that LEFM must be supplemented with information 
from smaller scales in order to make predictions about crack velocities, we should 
explain how this theory has been employed to make predictions about the limiting 
speeds of cracks. One can use LEFM to compute the energy travelling from far 
elastic fields in to the crack tip assuming that the crack is travelling at velocity v. For 
velocities v lower than the Rayleigh wave speed, this energy flux is always positive. 
For velocities above the Rayleigh wave speed, however, it is either imaginary or 
negative. There is one special case, of shear cracks travelling at y/2 times the shear 
wave speed, where the energy flux just vanishes, and this is the physically allowed 
value. Cracks in tension, by this logic, are forbidden to travel faster than the shear 
wave speed. 

Nevertheless, Buehler, Abraham, and Gao [59 ] observed such cracks in molec- 
ular dynamics simulations, and Petersan et al. 111051 observed them in experiments 
in rubber. They evade the apparent limitations of linear elastic theory because they 
correspond to a completely different scaling theory. 

In LEFM, the scaling parameter is 8c, the total displacement of the bound- 
aries at which the energy reaches the Griffith threshold and there is enough elastic 
potential energy available per unit length to create new surface. There is another 
critical extension, however, 8 C , which is the extension at which the solid would be 
stretched so much that the bonds between adjacent material points would give way 
completely. For most brittle materials, reaching 8 C requires stretching the solid to 
around 20% more than its original length. In rubber, it requires stretching the solid 
to 8 or 10 times its original length. 

If a crack tip avoids going unstable, then eventually when the boundary exten- 
sion 8 reaches a fraction of 8 C the crack speed becomes supersonic. The signature 
of the supersonic states lies in their behavior as the size of the system increases. 
For conventional subsonic cracks, cracks in two systems behave in the same way 
when £>g is the same. Thus if one strip is twice the height of another, cracks will 
have the same velocity in the two strips if the second is stretched by a factor of 
\/2 times as much as the first. For supersonic cracks, strain itself is the controlling 
variable, and energies rise to a new scale. Thus in the supersonic case cracks in the 
second strip will look like those in the first when it is stretched to a height 8 twice 
as much. Fig. [8]shows both both analytical calculations and experiments in rubber 
where these scaling behaviors are observed. The analytical calculations are carried 
out in exactly solvable lattice models II 1061 and the experiments are carried out 
in rubber. The calculations in the figure are for a simple case of anti-plane shear 
(Mode III), but they have been extended to fracture in tension as well 111071 . Rub- 
ber may seem an odd experimental setting for fracture mechanics, but the relation 
between force and displacement is nearly linear over wide range of displacements, 
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Supersonic cracks, scaled by extension 



Figure 8: (A) Four different views of Neo-Hookean crack velocities, showing that, 
depending upon how they are scaled and displayed, one focuses either upon con- 
ventional subsonic fractures or supersonic ruptures. In the limit of infinite system 
height N, the two different types of solutions are separated by an infinitely long 
plateau at the wave speed. The final panel simply re-plots the data of the third 
panel with a different range of velocities visible, emphasizing that if the macro- 
scopic limit is taken holding strain constant, it appears that all cracks are super- 
sonic. (B) Experiments in rubber at 85 C. The temperature needs to be elevated 
because at room temperature rubber undergoes a transition above extension ratios 
of 3.5 where its fracture energy increases enormously. The experiments are carried 
out in strip geometries with four different system heights. Velocities are scaled 
by the wave speed c = 21.9 m/s. The dimensionless extension 8/8g is estimated 
by the square root of the energy scaled by the Griffith energy of 5580 J/m 2 . The 
vertical extension ratio Xy, the ratio of the stretched height to the original height, is 
scaled by 10, which is a rough estimate of the critical extension 8 / 8 C obtained by 
comparing slopes in right and left panels. Experimental results courtesy of Hepeng 
Zhang, Johnathan Niemczura, and K. Ravi-Chandar. 
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and the strongest deviations from linear behavior are confined to a small region 
near the tip as required by small-scale yielding. 

The fundamental reason that cracks normally travel slower than all wave speeds 
is that they are only able to propagate by absorbing all the potential energy stored 
in the sample, all the way out to the boundaries. In the case of supersonic cracks, 
the solid has been stretched so much that energy sufficient to snap bonds is located 
with a fixed distance of the tip. Thus the crack speed is no longer limited by 
the time needed to transport energy to the tip from far away. Continuum elastic 
solutions for supersonic cracks show that their tip must now be a source of energy, 
a fact that has been used to argue that they cannot exist. In fact, enough energy 
is stored near the tip of a supersonic crack that it can break bonds with enough 
left over to travel outwards. A final observation is that in samples strained enough 
to support supersonic cracks, most of the energy in the sample is released through 
contraction of material after the crack has passed. This energy must be absorbed by 
some mechanism; in rubber, it is absorbed by bulk dissipation, which is particularly 
strong along a shock line that develops in the wake of the crack. 

4 Conclusions and Unsolved Problems 

We took a careful look at the underlying assumptions that describe our current view 
(LEFM) of how the simplest cracks evolve, and noted a few possible "cracks" in 
its underlying structure. 

• LEFM is a scale-free theory, but a fundamental understanding of the physics 
of fracture requires description at the scales where rupture processes occur. 

• Fracture processes are usually not explained by thermodynamic ground state 
arguments, and they are not even necessarily controlled by flow of energy to 
the crack tip. 

• The major unsolved problems of fracture concern the way the crack line de- 
cides how and where to move. Problems include the reasons cracks change 
direction while in motion, and conditions for crack fronts to remain essen- 
tially planar or to create complex ramified surfaces. 

The resolution of these issues may be intimately related to both our understanding 
of the proper "ground state" of this system and its stability. We have seen that there 
may not be a single solution to this problem, and different classes of physically vi- 
able solutions may exist. We surmise that a fundamental description of the stability 
of such solutions for "simple" cracks could be closely linked to how to properly 
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take into account the scales that define them. These scales may be dynamic en- 
tities, as we have seen, in sections 13.21 and 13.31 or static ones, as in sections 13.41 
and 13.51 They may simply determine the character of the ground state, or possess, 
themselves, intrinsic complex dynamics that may undermine ground-state stability. 
These issues define current active research directions. Perhaps their resolution (in 
the coming decade) will provide the basis of our next review. 
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